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p ^ We calculate a series of synthetic X-ray spectra from outflows originating 

O from the obscuring torus in active galactic nuclei (AGN). Such modeling includes 

2.5D hydrodynamical simulations of an X-ray excited torus wind, including the 
effects of X-ray heating, ionization, and radiation pressure. 3D radiation transfer 
calculations are performed in the 3D Sobolev approximation. Synthetic X-ray 
<^ line spectra and individual profiles of several strong lines are shown at different 

O inclination angles, observing times, and for different characteristics of the torus. 

^ Our calculations show that rich synthetic warm absorber spectra from 3D 

^ modeling are typically observed at a larger range of inclinations than was pre- 

2 viously inferred from simple analysis of the transmitted spectra. In general, our 

results are supportive of warm absorber models based on the hypothesis of an 
I— I "X-ray excited funnel flow" and are consistent with characteristics of such flows 

^ inferred from observations of warm absorbers from Seyfert 1 galaxies. 
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O 1. Introduction 

^ Among the most striking data taken with the grating instruments on the Chandra and 

^ XMM — Newton X-ray satellites are the rich line absorption and emission spectra from 

Seyfert galaxies and other active galactic nuclei (AGN). These spectra, referred to as 'warm 
absorbers', contain features from partially ionized ions of most astrophysically abundant 
elements, and they display systematic broadening and blueshifts corresponding to Doppler 
velocities ~ 10^ - 10^ km s 
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The blueshifts indicate an outflow, suggest mass fluxes high enough to be important to 
the AGN mass budget, and they may contain clues about the nature of the gas accreting 
onto the central black hole and its influence on the surrounding galaxy ( Murray et al. |[T995 ) . 

The X-ray and UV spectra are rich with spectral features and of high statistical quality, 
and appear to be directly associated with the orientation effects inferred from Seyfert galaxy 



unification (Antonucci & Miller 1985). That is, Seyfert 1 spectra show line absorption and 



scattering, while Seyfert 2s show line and bound-free continuum emission. The properties 
of the gas, i.e. the ionization, speed, and column density, are similar in the two cases; the 
primary difference is in the orientation relative to our line of sight. In Seyfert 1 objects we 
view the partially ionized gas primarily in transmission, and in Seyfert 2 objects we see it in 
reflection. 

Observed warm absorber outflow velocities, and virial arguments, suggest that they are 
launched at > 10^ gravitational radii, or ~ lO^^cm from a 10^ Mq black hole. This coincides 
with estimates for the location of the material, the 'obscuring molecular torus', which is 



responsible for blocking direct views of the central black hole in Seyfert 2 galaxies (Krolik & 



Begelma"ii||1986 ) . The structure of this obscuring torus, its origin and fate, are key questions 
in the study of AGN. Although the torus may represent part of a flow originating at smaller 



(Elitzur & Shlosman 2006) or larger (Proga 2007) radii, a convenient definition is that of 



cool (~ lOOOK) optically and geometrically thick gas in approximate rotational and virial 
equilibrium at ~1 pc (Krolik & Begelman 1988). Any gas within the central region of an 



AGN is likely to be heated and ionized by the strong UV-X-ray continuum from the black hole 
and inner accretion disk. The temperature of X-ray heated gas asymptotically approaches 
the 'Compton temperature'. Tic ~ lO^K, when the radiation flux is strong enough to exceed 
radiative cooling driven by atomic processes (in case of the Compton heated accretion disk 



winds see e.g. Begelman, McKee and Shields (1983)). However, adiabatic losses due to the 



wind motion may help to keep the temperature of the gas at lower levels (Dorodnitsyn et al. 



2008 b) 



Gas at the Compton temperature will not be bound in the local potential and will flow 
out in a wind (Krolik and Begelman 1986). The torus interior gas must be much cooler, 
and an X-ray illuminated torus will develop a heated intermediate region, corresponding 
to the transition between the two temperature extremes. The details of this transition are 
key to the understanding of warm absorbers, since neither the cool interior, nor the hot 
Compton flow can produce the spectra we observe. The structure of the transition depends 
on the density structure and, therefore, on the dynamics. Simple estimates, based on the 
assumption that the gas remains in ionization and thermal equilibrium with the radiation 
field, and that the flow geometry is spherical, predict that the transition should be sudden, so 
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there should be httle gas in the intermediate 'warm absorber' state. This is clearly in conflict 
with the observations, suggesting that time-dependent or geometrical effects are important: 
warm absorbers may represent a snapshot of this gas as it is transiently heated and flows 
out from the torus. 

We have constructed 2.5 dimensional (2D, axisymmetric with rotation) time dependent 
models for outflows driven by X-ray heating, which incorporate the torus geometry. Simple 
transfer solutions suggest these can make spectra similar to those observed. They produce 
column densities, mean ionization states, and velocities of the outflows which are comparable 
to what is inferred from observed X-ray spectra. These conclusions are backed up by de- 
tailed multi-dimensional numerical hydrodynamic calculations, incorporating the physics of 



radiative heating, cooling, photoionization, and radiation pressure (Dorodnitsyn et al. 2008 
afb ). A challenge for these models is that, since the warm absorber gas we observe is in a 
transient state, it occupies only a fraction of the volume in the torus opening. This predicts 
detectable warm absorbers from a smaller fraction of AGN than is observed, since for many 
inclination angles the line of sight to the X-ray continuum source misses the partially ionized 
gas. The observed fraction is ~ 50%, and perhaps greater (Reynolds, 1997), while models 
predict approximately half this number. Furthermore, the models calculated so far do not 
reproduce some of the features observed from X-ray spectra, including an apparent bimodal 
distribution of ionization in many objects, and the presence of X-ray absorption by dust or 
neutral material in some spectra. 

Some of the conclusions of our past work are limited by the simple 1-dimensional transfer 
we employed to calculate the emergent spectra. Although these were useful as a preliminary 
step of predicting line spectrum, the ID radiation transfer calculations are not capable of 
producing emission lines and of treating effects like fllling of absorption lines with emission 
lines. In addition, all effects associated with non-radial motion of the fluid were neglected. 

In this paper we produce full 3D radiation transfer calculations of X-ray line spectra 



from such flows using the information provided by our hydrodynamical models (Dorodnitsyn 
et al.||2008 alb) ). These simulat ions produce spatial distributions of the projected velocities, 
and densities. Using this information as an input to the 3D radiation transfer in X-ray lines, 
we perform 3D radiation transfer in the Sobolev approximation. We calculate spectra in 
approximately 10^ lines and also high resolution profiles of some of the strongest lines. 
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2. Methods 

To calculate an X-ray line spectrum from the moving gas which was stripped off the 
AGN torus via X-ray induced evaporation, we combine together the following ingredients: 

i) Taking into account heating and cooling processes we calculate a set of time-dependent 
2.5D hydrodynamical models of winds evaporated from the AGN torus, ii) Based on the 3D 
generalization of the escape probability method, we construct a numerical code which is able 
to calculate line spectra from a 3D flow, iii) Making use of the XSTAR code we calculate 
tables of X-ray opacities for different energies and ionization parameters, iv) Using opacity 
tables and distributions of density, and velocity from hydrodynamical models, we calculate 



the emergent spectrum. Step i) has been described in Dorodnitsyn et al. (2008 a). In the 



following we briefly review methods and describe results of steps ii) through iv) 



2.1. Hydrodynamical modeling 



Our hydrodynamic and radiation framework is described in detail Dorodnitsyn et al. 
( |2008 a||b| . Here we briefly review the methods and results. 



The primary physical assumption of our model is that that the inner face of the rota- 
tionally supported torus is exposed to radiation heating by X-rays generated by accretion 
near the supermassive black hole. 



The equations of hydrodynamics are calculated using the ZEUS-2D code (Stone & 



Norman 1992), modified to take into account various processes of radiation- matter inter- 



action; this modeling includes time-dependent 2.5D hydrodynamical simulations (that is 
time-dependent 2D simulations with rotation) of the flow. Interaction of the radiation with 
matter is incorporated into the hydrodynamical code taking into account radiation heating 
and cooling in the energy equation and pressure of the radiation in UV spectral lines and 
in continuum in the momentum equation. The rates of Compton and photo-ionization heat- 
ing and Compton, radiative recombination, bremsstrahlung and line cooling are calculated 



using approximate formulae, which are modified from those of Blondin (1994) (Dorodnitsyn 



et al. 2008 b). These modified formulae adopt newly available atomic data and heating and 



cooling rates obtained from the XSTAR code (Kallman & Bautista 2001). 



We assume that a certain fraction of a total black hole (BH) accretion luminosity, is 
emitted as X-rays which interact with the inner face of the rotationally supported torus 
and heat it. Thus the cold gas of the torus is heated to almost the Compton temperature 
(T ~ 10''), fills up a throat of the torus constituting an outflow. This hypothesis was tested 
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both qualitatively and quantitatively by evolving such models in time, and exploring the 
parameter space, including the accretion efficiency, and geometrical and physical properties 
of the torus. 



2.2. Radiation transfer in spectral lines 



In the papers by porodnitsyn et al.| (2008 a|[b| pure absorption spectra were calculated 



from different models assuming different inclinations. Transfer was calculated using simple 
ID attenuation of the radiation due to line and bound-free absorption, assuming pure radial 
streaming of the radiation and also approximating the emitting core as a point source. 
The obvious limitations of this approach motivates a much more realistic radiation transfer 
treatment which would take into account spatial gradients of the velocity in a 3D flow. In 
this work we implement such an approach, in a 3D radiation transfer calculations in Sobolev 
approximation for the line radiation transfer. 

The supersonic motion of gas of the AGN wind provides justification for use of the 
Sobolev approximation while calculating radiation transfer in spectral lines. The problem of 
radiation transfer in spectral lines in the Sobolev approximation in a spherically-symmetrical 



flow is addressed in numerous papers (for references see e.g. Rybicki & Hummer ( 1978 , 1983 ); 



Dorodnitsyn (2009) ). The primary assumption of the Sobolev approximation is that due to 
velocity gradients in the flow the Doppler shifting of frequency moves a photon quickly out of 
the resonance soon after it interacts with the flow. Thus the radiation transfer in the Sobolev 
approximation is treated as an entirely local problem, assuming the characteristics of the 
absorbing plasma are constant across the narrow region of interaction. Our hydrodynamical 
model is axially symmetric, but this symmetry is lost if viewed from an inclination angle, so 
to calculate the spectrum we must adopt a 3D Sobolev approximation. Even if the poloidal 
velocity along the streamline is monotonically increasing, the projected velocity may not 
behave in such a way, rather having a maximum or several maxima (depending on the real 
structure of the 2D flow). This poses the problem of multiple resonances which we discuss 
in detail below. 

The geometry of the problem is shown in FigjlJ The Cartesian coordinate system xyz 
{y axis is not shown) is located at O. The observer is situated at inflnity looking at the torus 
from an inclination, ^o- Since the hydrodynamical model is axially-symmetric, we put the 
observer in the xz plane. The plane of the sky, PS, is transported to O, providing another 
Cartesian coordinate system x'y'z', [y' axis is not shown) which is the xyz system, rotated 
by the angle around the y axis. Notice also that y' = y. Our assumptions imply (see 
further in the text) that a photon after being scattered in a spectral line, for example, in the 



point R, propagates along the straight hne RB in the direction of the unit vector n. This 
ray, RB intersects the plane PS at the point A, being at the angle 0o from the OZ' axis and 
at a distance p (impact parameter) from O. Summing all such possible rays, i.e. integrating 
over the plane, PS, the radiation flux which is observed at infinity is obtained. 



Fig. 1. — Illustration of the geometry and coordinate systems implied by the calculations 
(see text for details). Notice, that the y-axis (not shown) points into the plane of the picture 
to form a right-handed xyz coordinate system, and y' = y. Not to scale. 

The total radiation flux which is registered by the observer at infinity at the frequency, 
equals to 



Jo Jo 

where Fc is the flux in the continuum near the line frequency. 

In our treatment of emission only resonant lines are considered. Although we do consider 
non-resonant absorption in the approximate treatment of the attenuation (section 2.3). After 
being emitted by the core, or, alternatively, being conservatively scattered in a line, a photon 
remains at constant frequency, z/°° and direction of propagation, unless its frequency, Ucom 
in the frame co-moving with the gas appears in a resonance with a line transition: 



Z 





Toms 




(1) 
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where z// is the frequency of the hne, and !S.vd is the Doppler width of the hne. For such 
a z/°°, all such possible resonant points constitute resonant surfaces, which could have an 
extremely complicated shape in the case of a real flow. After encountering the resonance 
surface the intensity J changes discontinuously: 



where Si = S{ai) is the source function and Ti^i = r(cri) is the optical depth in a spectral 
line. The intensity changes discontinuously, being a) attenuated by a e"'^'''*^"''-* factor, and b) 
reinforced by the contribution of the locally scattered radiation field, Jioc,i = >S'i(l — e""^' ') ~ 
l/r^(l — e""^''^) (the last equality holds true for a point source and pure line scattering). The 
latter contribution results from scattering of external photons at a resonant point. Let s be 
the distance measured along the photon's trajectory, then r; in a direction n in a moving 
gas depends on the gradient of the velocity in this direction, i.e. on dv/ds and can be cast 
in the form: 



P e KesVth 



(4) 



\dv/ds\ ' 

es, Hi is the line opacity, Kes is the Thomson scattering opacity, p is the 



where e = ki/k. 
density, and Vth is the thermal velocity. 

In the presence of multiple resonances, distant points within the flow are causally con- 
nected, i.e. they can exchange photons ranging within the line frequency width. In terms of 
surfaces of equal frequencies, such surfaces become multi branched, allowing for a photon of 
a frequency to interact with the moving matter at different points along a single line of 
sight. In the presence of multiple resonance surfaces, we assume that the interaction between 
crj_i and cTj takes place only in the forward direction (in the direction towards the observer). 
All other possible interactions are neglected. Then the expression (jsj) should be just summed 
over the available resonance surfaces. Altogether, these assumptions are known as the dis- 



connected approximation (Grachev & Grinin 1975; Marti & Noerdlinger 1977 Rybicki & 



Hummer 1978, 1983): 



N 



r°°(z/°°l 



jins.O g- Y.'i' ni 



+ E (1 - ^^"0 



- nj 



(5) 



where N is the total number of resonances and / 



ins,0 



-^src e"'^"^'-'""^ is the initial intensity 
of the radiation at the point Tq, which is the farthest point along the ray from the observer 
in the computational domain, Jsrc is the intensity emitted by the source (core), and for 
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simplicity, we omit the dependence of / on p, and 6^0; t^s is an electron scattering optical 
depth: 



Tes{r) = / Kespdr. (6) 



The source function is calculated with the method of escape probabilities (e.g. Castor 



(1970)). The probability of a photon to escape in a direction n in a moving gas, depends on 



the gradient of the velocity in this direction, n, i.e. on dv/ds 

dv 



J]ei,jn,nj, (7) 



where is the rate of strain tensor: 



e.J = UI^ + l^'l, (8) 



2 \drj dvi J ' 



Components of (\8v in spherical polar coordinates, are given, for example, by Batchelor 



(1967). 



A probability, of a photon to escape in a direction n is given by 

Pn = {l-e~-^-^)/e-<-\ (9) 

Integrating this relation over solid angles Q{n) (i.e. either all over the sky, or over solid angles 
subtended by the source, flsrc) results, respectively, in escape and penetration probabilities 



(Rybicki & Hummer 1983) 



Pesc = ^ I Pndn. (10) 



An 



/^pen= 73 / (^ndn. (11) 



src 



4tt 

In the case of pure line scattering, the source function is cast in the form 



Sir) = (12) 

Pesc 



where J is the mean intensity at the point r: 
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J(r) = ^/" f3X,{r)dQ, (13) 

and /srclr) is the intensity of the continuum near Pi, possibly attenuated by the electron 
scattering prior to resonant point, r: /src(r) = -^srcC""^". It is assumed that the attenuation 
is dominated by Thomson scattering k^s = 0.2(1 + Xh) — 0.4cm^g~^, where Xh is the mass 
fraction of hydrogen, and the factor e""^*" , accounts approximately for the attenuation of the 
radiation flux on the way from the source to the resonant point. 

Pre-calculated hydrodynamical models provide the distributions of density, p(r) and 
velocity, v(r), needed to calculate Si and th. With this knowledge in hand, substituting 
I°°[u°°) from ^ to ([T]), we calculate the spectrum. Notice that in case of spherically 
symmetrical distributions of the velocity, density, etc., the above formalism gives the well 
known P-Cygni profiles. 



We have rigorously tested our calculations against the atlas of P-Cygni profiles of Castor 



& Lamers (1979). Results of such testing are shown in Dorodnitsyn (2009), Figure 9. 



2.3. Treatment of lines and continua 



First we point out that in this work we consider transitions in all lines as resonant. 
That is, we don't differentiate between true resonant event and lines whose upper levels are 
depopulated by photoionization, auto-ionization, and thermalization processes. This is an 
approximation which has only a small effect on our final results, since the strongest features 
in our spectra are true resonance lines. Using the XSTAR code we calculate two opacity 
tables: table 1 contains only line opacities, and table 2, only continua opacity. Both tables 
contain cross sections as functions of ionization parameter (see further in the text) and 
energy. As mentioned in section 2.2, to integrate flux, F, from ([T]) we sum all rays with all 
impact parameters, p at all < < 27r (see Figure [T|. Starting from the boundary of the 
computational domain (which is at r = 50pc), we calculate intensity along the ray RB. At 
each point we calculate the local Thomson optical depth between this point and the nucleus 
(located at r=0.05pc), then calculate the attenuated ionization parameter, and then refer 
to table 1 to check the resonant condition, (|2|, and if succeeded calculate the new value of 
intensity, /, using ([3]). Between resonant points we attenuate the radiation field, making use 
of pure continuum opacities, from table 2. We do not include any turbulent broadening in 
any of our calculations, rather just thermal broadening plus the effects of the bulk flow via 
the Sobolev optical depth. 
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2.4. Assumptions and input parameters 



Our hydrodynamical modeling is based on the solution of the 2D, time-dependent hy- 
drodynamical equations for rotating flow exposed to external X-ray heating. Radiation 
pressure is taken into account both in continuum and due to UV pressure in spectral lines 



in which case it is calculated from the generalization of the Castor et al. (1975) formalism 
dPorodnitsyn et al.|[2008 aP . 



The assumed photo-ionization equilibrium of the wind implies the outflowing gas can 



be parameterized in terms of the ratio of radiation energy density to baryon density (Tarter 



et al. 1969). One popular form of this parameter is 



e = 47rFjn, (14) 

where Fx = Lxe~'^/(47rr^) is the local X-ray flux, Lx is the X-ray luminosity of the nucleus, 
and and n is the number density. We assume that the nucleus or the "core" at radius 0.05 
pc, emits a continuum power law spectrum with energy index a = 1. 

The ionization parameter ^ can be estimated in the following approximate way: ^ ~ 
4 ■ 10^ ■ /x r Me/(A'"23 Tpc), where A^23 is the column density in 10^^ cm~^, /x is a fraction of 
the total accretion luminosity Lbh available in X-rays, F is a fraction of the total Eddington 
luminosity Ledd = 1-25 • 10^^ Mg, where Mg is the mass of a BH in IO^Mq. The radiation 
pressure in lines is calculated assuming the relative fraction of UV radiation, /uv =0.5. 
Heating and cooling rates account for Compton and photo-ionization heating and Compton, 
radiative recombination, bremsstrahlung and line cooling, and are incorporated into the 
equations of hydrodynamics. 

After non-dimensionalisation, our problem includes several characteristic scales: the 
radius is measured in terms of Rq, which is the distance of the initial torus density max- 
imum from the BH; time in terms of = R^J"^ / y/GM ~ 4.5 ■ 10^^ rpc^ Mg"^''^ (s), where 
Tpc is the distance in parsecs; the characteristic velocity reads Vq = a/ GM/ Rq ~ 6.6 ■ 



Vpc^'^ (cms ^) . From now on we measure time, t in units of to- 
Initial conditions at t = are provided by the equilibrium toroidal configuration in 
which a polytropic, rotating gas is in equilibrium with the force of gravity and radiation 



force. Thus, the initial state is a rotationally supported torus of Papaloizou fc Pringle 



(1984), which is corrected for the radiation pressure. The parameter that determines the 
relative distortion of such a torus is d = {w~ + w^)/2, where w is the cylindrical radius, 
and the inner and outer edges of the torus are located at w~ and w"*", respectively. 
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3. Results 

The major goal of this paper is to calculate emission and absorption spectra predicted 
by hydrodynamical simulations. The maximum density of the initial torus, rimax; or, al- 
ternatively, its initial Compton optical depth = r(6' = 90°) is used as a parameter to 
distinguish between different models; nmax roughly scales with the column density of the 
absorbing gas during the evolution of the torus. 

Our modeling of the spectra use a set of pre-calculated hydrodynamical models, includ- 
ing combinations for: = 1.3 (models Ai), and = 40 (models Bi); equivalently, these 
correspond to initial tori, having nmax = 10'' cm^^ or nmax = lO^cm"^. Other parameters, 
determining the torus arc-Ro and the distortion parameter, d. Models Ai and Bi have d = 2.5 
Additionally there are two models with d = 5 (models Ci) which represent an extreme case 



of extended and rarified torus. The full list of models is given in Table 1 of Dorodnitsyn et 



ah (2008 b). 



Spectra calculated from hydrodynamics contain too many features for easy interpreta- 
tion, and to make a comparison with observations feasible we convolve them with a Gaus- 
sian which has FWHM(E)= 10~^-E'; such convolved spectra are most closely related to what 
would be observed by an instrument such the Chandra grating. However, contrary to such 
a real observer, we are also aware of un-convolved spectra, and we frequently use these to 
compare and interpret the results of our analysis and to establish possible pitfalls. Further, 
we show that the interpretation of individual lines from the convolved spectra may contain 
significant caveats leading to a misinterpretation of the positions of some lines and of their 
relative blue-, and redshifts. Our spectra contain absorption as well as emission features 
and for convenience we denote the wavelength (or, equivalently the energy) of an absorption 
line as A and of an emission line as A; the line centroid as Aq; the blue- ward edge as A^, 
and the red- ward as A~. In this paper when speaking about individual lines we adopt the 
following convention: velocities corresponding to blueshifts are positive and to redshifts are 
respectively negative. 



From hydrodynamical modeling ( Dorodnitsyn et al.|2008 a ), it is established that by the 



time t ~ 3 the torus is losing mass via a quasi-stationary evaporative wind. Again, times are 
measured in terms of t^. Before this time, say at t ~ 1 the influence of the initial switching 
on of the X-ray heating is felt by the wind -|- torus system, and at much later time, say at 
t ~ 6, the torus becomes significantly influenced by the process of evaporation. Thus for our 
investigation we choose evolution stages at which the torus possesses a well developed wind 
but still is not drastically influenced by it. The results of our modeling for models B^, B^, 



and ^3 of Dorodnitsyn et al. (2008 b) at t=4 are shown in Figures 2 - 5 
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Results for the Bq model 
Results for this model at t=4 are shown in Figure [2j The initial state has Rq = 1 pc, 



40 n 



10" cm ^, and d = 2.5. 



Characteristics of the individual lines are consistent with hydrodynamical models. At 
6 ~ 35° mostly emission lines are seen atop of the continuum spectrum and at ^ = 51°, emis- 
sion lines are observed atop a significantly absorbed continuum. The position of a blueshifted 
O VIIIKa, 18.97A(653.5 eV) line remains approximately constant at all inclinations. For ex- 
ample, at 6' = 45° the velocity of the gas is —725 < f < 880 km s~^; the line centroid is 
Ao = 18.96A (653.7) eV, corresponding to f ~ 55 km s^^, which is comparable with the 
escape velocity at 1 pc. 



Rest Wavelength (Angstroms) 



40 23.84 16.98 13.19 10.69 9.05 7.85 6.93 6.2 




^.UaA.iI. 



0.8 1 1.2 

Rest Energy (keV) 



Fig. 2. — Model Bq (r^ = 40, Rq = Ipc) X-ray spectra, observed at t = 4 as a function of 
angle. 

At ^ = 47° a strong absorption line is observed at Aq = 16. 8A (738.3) eV in the convolved 
spectrum. If this line is due to OVII, at 17. 2A (720.84 eV ) this transition indicates an 
7198 km s~^ outflow. The escape velocity at distances which are plausible for the formation 
of such absorption features (~ 1 pc) is significantly smaller, ^ 50 — 60kms~^, implying such 
an interpretation is unlikely. In fact, looking at the un-convolved, high resolution spectrum 
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Distance along the line of sight (pc) 



Fig. 3. — Model Bq, observed at time t = 4 and at inclination 6q = 45° as a function of an 
impact parameter, p; the angular coordinate in the plane of the sky, 0o = 45°. Left panel: 
column density; right panel: ionization parameter. Legend: each curve is calculated for the 
impact parameter, p, measured in parsecs. 
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reveals a complicated picture in this energy range ~ 17.2 — 16.7 A (720 — 740 eV), consisting 
of many absorption lines which are blended together in the convolved spectrum. From the 
un-convolved spectrum we attribute a Aq = 15. 2A (722.1 eV) feature in the un-convolved 
spectrum to O VII, 17.2A (720.84 eV) line, blueshifted with v ~ 524kms~^. The position 
of the OVII, 16. SA (739.3 eV) edge, clearly seen in the un-convolved spectra indicates v ~ 
700 km s"^. This is indicative of a strong variation of the ionization parameter, ^ along 
the line of sight in the spectrum forming region (compare to Figure [3] or to Figure 4 from 



Dorodnitsyn et al. (2008 b)) 



The attenuation of the radiation field at high inclination is consistent with increasing 
column density, N{9) at high 9. The radial column density rises from N22 — 1 at r = 0.93 pc 
to N22 = 3.3 at r = 2.9 pc up to N22 = 75 at r = 3.9 pc, where N22 is the column density in 
terms of 10^^ cm~^. 

The line centroid, Aq is determined by two regions, where the absorption is most fa- 
vorable: starting, approximately from r = 0.8 pc, where — ^ 115kms~^ and den- 
sity is n ~ 4 ■ lO^cm"^; at r = 1.8 pc the density, n slightly drops, but the velocity 
rises: Vr — 300kms~^; at r = 3.56 pc the density has a maximum on the line-of-sight: 
n ~ 2.7 ■ lO^cm"^, but the velocity is only Vr — 26kms~^. The ionization parameter, ^ has 
a plateau, where ^ ~ 700, then it behaves like inverted number density, i.e. drops to zero at 
r ~ 3.5 pc, where n in turn has a maximum. 

An apparent emission hump at ~ 14.9 — 17.2A, (720 — 830 eV), 9 = 47° is due to 
H-, and He-like oxygen and moderately ionizated iron, FeXVII — FeXXIV. It has an onset 
near the OVII edge at 16.8A (739 eV) and is bounded by the OVIII edge at 14.3A (870 
eV) from the high energy side. At higher inclinations this emission excess is substituted 
by the absorption trough seen roughly at the same frequencies. As shown in Figure 4 of 



Dorodnitsyn et al. (2008 b) it is approximately this inclination where there is a significant 
rise in column density and the line of sight goes though the dense and low ionization gas. 

From the strong absorption line of MgXI at 6* = 47° at Aq = 9.2A (1353 eV), we infer 

V ~ 221 kms~^; from the high resolution spectra we get —665 < v < 887 km s~^. At smaller 
inclinations this line is not prominent; the MgXII line 8.4A (1472 eV) is most prominent at 
9 = 47° giving Aq = 8.4A (1473 eV), —407 < f < 611 km s~^, and also from the line centroid: 

V ~ 203 kms^^ 

The emission spectrum contains many lines. For example, emission is seen close to 
oxygen edge, e.g. OVII, indicating v ~ 162kms~^. At slightly lower energies many lines are 
blended together, producing a feature ~ 0.3A (10 eV) wide. Among other features, most 
prominent are those from MgXI, MgXII, NeX. Blueshifts of these lines are consistent with 
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a wind which is moving with velocity, 50 < f < 600kms~^. Positions of emission hnes, i.e. 
A~, Aq, and A"*" are consistent with that inferred directly from the hydrodynamical models 
(i.e. for example with those given in Figure 4 of Dorodnitsyn et al. (2008 b)). 

At = 10°, the maximum projected velocity is determined by the maximum v.,. ^ 
664kms~^ at r ^ 3.3 pc. The toroidal velocity, peaks at r ~ 0.37 pc at a value ^ 
500 km s"'^. The column density smoothly rises from N22 = 0.15 at r ~ 1 pc to N22 — 1-8 at 
r = 20 pc. 

Comparison of the spectra given in Figure |2] with those given in Figure 10 of Dorod-| 
nitsyn et al. (2008 b) (which shows the pure transmitted spectra for the same model Bq) 



demonstrates significant differences. That is, the full 3d transfer models allow for more trans- 
mitted flux at a given viewing angle. These differences are partially due to the higher energy 
resolution of the opacity tables adopted in the current simulations and partially because of 
the qualitative difference of the 3D approach (presence of significant emission component, 
finite size of the core). 



Obscuration of the emission lines 

No significant negative velocities are detected in emission lines. This is the result of 
a severe obscuration at high impact parameters, p. Let L be the distance between points 
A, and R, Figure [l] In Figure [3] is shown the behavior of the column density, N{L,9q,(I)q) 
and the ionization parameter, ^(L,6'o,0o) with distance L. Figure [3] shows these quantities 
at = 45°, 00 = 45°. The overall behavior indicates a significant absorption at high p. 
What we see is the line-of-sight penetrating more cold gas as it moves from p = 1 to p = 5. 
Then, at larger p = 10 there is again less gas (this is the end of a torus), and reducing to 
smaller values at p = 30. Note that if there is only resonant scattering in lines present then 
significant emission can be obtained by accumulating photons form a large volume, which 
has a small Thomson depth. Our result shows that it is difficult to see extended regions 
behind the torus, and thus no significant emission lines with negative velocities are seen in 
the simulations. 



Results for the model 

This model represents a wind which is launched closer to the BH, at Rq = 0.5 pc. The 
overall spectra are similar to that of Bq model although the larger escape velocity, Vesc{Ro) 
implies larger velocities inferred from individual lines. 
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The prominent O VIIIKa line is blueshifted: Aq = 19A (653.8 eV), v ^ 100 km at 
9 = 45°. At higher inchnations, 6 = 47°, it is stronger but remains with the same blueshift; 
at 6* = 51°, Ao — 19A. Such blueshift is consistent with the gas moving at Vcsd^pc). At 
higher energies, among strong absorption lines are those of NeX, f ~ 500 km s-^ MgXI and 
MgXII, indicating \v\ < 1000 km s~^. 
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Fig. 4. — Model B3 (r^ = 40, Rq = 0.5pc) X-ray spectra, observed at t = 4 as a function 
of angle. Note a different scale in the y in the plot at the bottom; at ^ > 55° the X-ray flux 
is strongly suppressed by the Compton thick cold matter. 

At higher energies. He-like Si Ka line is prominent (Figure |4|: Aq = 6.6A (1866 eV), 
t>(Ao) — 160kms~^, v{X^) ~ 321kms~^, f(A^) ~ — 161kms~^ At low/high inclinations, 
when the emission is seen on the background of either absorbed or unabsorbed continuum, 
an excess emission bump is seen around 16. 8A (740 eV); many lines which are close to O VII, 
and O VIII are blended together contributing to this emission feature. Observing the wind 
from inclinations which are "just right" (in our case ~ 45°) allows us to see the absorption 
from the gas which is just on the edge of the torus; this is the case when absorption lines 
may show less blueshift than emission. Some emission comes from the gas which is at lower 
6 which is generally moving faster than the wind at high inclinations. From Si XIII we have: 
—800 < V < 805kms~^ from 6 = 35° emission, and —650 < v < 600kms~^ from 6 = 45° 
absorption. We conclude that results for model are consistent with the general situation 
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in which the gas is located two times closer to the BH, and the familiar funnel mechanism 
also helps it to escape with larger velocities: 50 < f < 1000 km s~^. 

The structure of the hydrodynamical flow of this model is similar to that previously 
discussed: the radial column density at 6' = 45° rises from N22 — 0.5 at r = 0.93 pc to 
N22 = 1-7 at r = 2.9 pc up N22 = 38 at r = 3.9 pc; remarkably, the gain in column density 
occurs at the same radii as in the Bq model. The line centroid, Aq as in the Bq model, is 
determined by two regions of the most favorable absorption: one starting from r ~ 0.8 pc, 
where Vr — ^ 162kms~^ and density is n ~ 3 ■ lO^cm"^ and the other at r ~ 3.56 pc 
where the density has a maximum value on this line-of-sight: n ~ 2.74 ■ lO^cm"^, and the 
velocity Vr — 37kms~^. Between these regions at r = 1.8 pc the density decreases but the 
velocity increases: Vr = 416kms~^, as expected from mass flux conservation. As in the Bq 
case, the ionization parameter, ^ has a plateau, where ^ ~ 2000. Then it behaves just like 
inverted number density, i.e. to very small values at r ~ 3.5 pc. 

At high inclination the radiation of the nucleus is highly absorbed by cold gas of the 
torus. There are a few emission lines seen on the background of such attenuated continuum. 
Those include Aq ~ 13. 5A (921.4) eV representing a NelX line, possibly blended with mod- 
erately ionized FeXVII — FeXXIV lines. From the broad NeX line, we obtain A+ ~ 12. lA 
(1024 eV), A~ ^ I2.2A (1020 eV), indicating \v\ < 587kms-^ From MgXI line, we obtain 
A+ ~ 9.I5A (1354 eV), A" ~ 9.17A (1352 eV), indicating \v\ < 443 km s'^ Maximum 
velocities of emission lines are reminiscent of the high velocity, low inclination outflow. At 
low inclination, 6 ~ 10°, the hydrodynamic model gives Vr — 935kms~^ at r ~ 3.2 pc and 
toroidal velocity, ^ 705kms~^ at r ~ 0.37 pc; the column density smoothly rises from 
N22 = 0.1 at r ~ Ipc to N22 ~ 1 at r = 20 pc. Overall attenuation of emission lines are 
consistent with the picture inferred for the Bq model attenuation (c.f. Figure [3]). 



Results for the model 

Figure [s] shows convolved emission-absorption spectra for the model A3 (rj^ = 1.3, 
Rq = 0.5, r = 0.5, d = 2.5) at time, to = 4 at different inclinations. A pure emission 
line spectrum is showing up at < 47°. Strong lines from OVII, OVIII, MgXI, MgXII, 
Si XII, are blended with many other lines, in particular with those from low ionization 
FeXVII — FeXXIV. Inferred velocities are consistent with those known directly from the 
hydrodynamic model. 

Observing emission lines at low 9 we see directly through the torus throat, thus having 
a chance to see the gas which recedes from us on the other side. Thus, at = 35°, from 
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Fig. 5. — Model A3 (r^ = 1.3, Rq — 0.5pc) X-ray spectra, observed at i = 4 as a function of 
angle. At ^ > 55° the X-ray flux is strongly suppressed by the Compton thick cold matter. 
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O VIII Ka line we infer A+ = 18.9A (655.5 eV) and A" = 19A (652 eV) indicating -725 < 

V < 881kms~^. Many emission lines at higher energies have symmetrical profiles: from 
MgXI line we infer A+ = 9.2A (1355 eV) and A" = 9.2A (1349 eV), indicating -665 < 

V < 665kms~^; at higher energies, from Si XIII line we obtain A"*" = 6.66A (1865 eV) and 
A~ = 6.64A (1861 eV), indicating —643 < v < 643kms~^. Observed from an inclination 

— 45° we look almost in the torus throat, just touching its inner evaporative skin, and 
these emission lines are formed by the gas within such a cone, i.e. at < ^ < ^o- We note 
the likely role of the accretion disk which in real AGN is opaque to the external radiation 
and can screen the receding gas. Obviously, no accretion disk in our simulations is explicitly 
taken into account so we are not able to take this effect into account. 

The lower density of models Ai in comparison with models Bi, allow us to see more 
interesting profiles of individual lines. At the critical inclination, 6 ~ 45° the spectrum 
resembles that obtained from B3 or Bq models. However, in this case the profiles of many 
lines have a shape of P-Cygni and many other have emission-absorption-emission features, 
with absorption being considerably stronger than emission. One example of such an M- 
shaped profile is shown in the next section (Figure [9]). The red- or blueshifted wings seen 
in the emission are due to a combination of rotation and the P-Cygni mechanism (in the 
red- ward part of the profile) and due to rotation (in the blue- ward part of the line profile). 
The presence of M or U shape profiles is clearly marking the problem: the importance of 
rotation and fast gas moving at low r, i.e. v^, and vg. Note that the appearance of such 
profiles from a rotating ring or disk, in the Sobolev approximation, has been emphasized by 



Rybicki & Hummer (1983). We will discuss these effects in some detail further in the text. 



Figure [6j shows the distribution of column density, N(L, 9q, 0o) and ionization parameter, 
^(L, 6*0, 0o) along the line of sight. At p ~ 1, L ~ 0, the edge of the cooler slow wind or torus 
can be seen in the increase of N. Along the line-off-sight the ionization parameter is in the 
range 1 < ^ < 10^. 



4. Discussion 

Analysis of different models shows that, roughly speaking, in the middle of the torus 
evolution the spectrum which is most rich with spectral features is observed at inclinations 
near 6 = 45°. This holds approximately true for both types of hydrodynamical models, Ai 
and Bi. At lower inclinations gas is mostly observed in emission lines while at higher incli- 
nations, the spectrum becomes significantly absorbed by the cold gas in the torus. However 
some of the strong emission lines are seen atop of such absorbed continuum. 
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Fig. 6. — Model ^3, observed at time t = 4 and at inclination ^0 = 45° as a function of 
impact parameter, p; the angular coordinate in the plane of the sky, 0o = 45°. Left panel: 
column density; right panel: ionization parameter. Legend: each curve is calculated for the 
impact parameter, p, measured in parsecs. 
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The positions of the hne centroid energies are generally consistent with what is expected 
from the torus X-ray exited wind paradigm, i.e. the inferred velocities are several times 
at 1 pc. There is also found an influence of the gas that was evaporated off the torus and 
then accreted deep into the potential well; larger covering factor and high velocity allows it 
to block enough continuum and produce high velocity broadening of some lines. These high 
velocity signatures come from distances as small as 0.1 pc. Another interesting effect is a 
blending of multiple absorption and emission lines which are seen in the un-convolved, high 
resolution spectrum but form a single absorption line in the convolved spectrum. Comparing 
high resolution and convolved spectra, we show that interpretation of such lines may in some 
cases lead to significant overestimation of the inferred gas velocity. 

Many strong lines are observed blueshifted in absorption and at almost zero blueshift 
or very small in emission: f(Ao) < 50kms~^; For example, MgXII line 8.4A (1472 eV) 
in absorption indicates v ~ 200kms~^ and having in emission almost zero shift, Aq — 
Ao(MgXlI) for the model B^. However from the high resolution spectrum we see the wings 
of this line indicate v ~ 600kms~^ (for the model) 

Blending of lines is important, for example there is a spectral feature which is observed in 
emission at ~ 50°, and in absorption at 6* ~ 45° in B-^ and Bq models. Comparing with the 
un-convolved spectrum shows that this is the result of blending of many lines, specifically the 
H- and He-like Ne, and lines from low ionization Fe occurring in the 13.4-13.5A (920.2-925.7 
eV) range. After convolution it produces a spectral line at 13. 5A (921.9 eV). 



4.0.1. M shape line profiles 



Individual line profiles obtained from some of our hydrodynamical models display multi- 
component features. That is, such a profile has two emission components separated by an 
absorption trough. Most clearly such profiles are seen in models Ai i.e. in those models which 
tend to have lower density and a wider throat of the torus. U shape profiles are well known 
in the case of fluorescent line emission from accretion disk (see e.g. the nonrelativistic limit 



of 


Cunningham 


(1975 


) or 


Gerbal & Pelat 


(1981)). 


Rybicki & Hummer 


( 


1983 


calculated a profile of a 



Contrary to the familiar reflection fluorescent line, in their case the line is formed due to 
resonant scattering within the disk with the resultant formation of an M - shaped line profile. 
Roughly speaking, motion of the gas distorts the M-shape profile towards a P-Cygni shape, 
i.e. contributing to the red-ward emission hump. Altogether there appears to be an M-shape 
profile with stronger red- ward emission (Figure [9]). 
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Fig. 7. — Model Aq, observed at different evolution times at inclination 60 — 47°. 
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4-0.2. Evolution in time 
The time evolution of model Aq is shown in Figure [7) 

The O VIII La is seen in absorption from t ~ 2 to t ~ 6. Starting from t = 2.5 to 
t ~ 6 from the high resolution spectra, we have: f(Ao) — 45kms~^, f(A~) ~ — 642kms~^, 
and f(A^) ~ 872kms~^. Such characteristics remain almost unchanged until t ~ 6, at 
which the weak M-shape profile appears: Aq remains approximately the same, while v{\~) ~ 
-550kms"\ f(A") ~ 640 km s"^ and w(A+) ~ 642 km s"^ t'(A+) ~ 870 km s"^ ( A=^ is the 
energy of the most red-, or blue- ward wing of the line, A^ is related to the different emission 
parts of the single M-shaped line). 

Quite different behavior is found for other lines, such as Si XIII: the line is in absorption 
at times t < 5 showing weak P-Cygni profiles at t = 5.5: f(A~) ^ — 643kms~^, f(Ao) = 
160kms~^, and f(A"^) = 321kms~^. 

During the time evolution the torus is heated and loses mass via the wind. Lowest 
density parts of it suffer most, and the torus is additionally influenced by the combined 
pressure of its own wind and radiation pressure. Altogether they work to squeeze the torus 



towards the equatorial plane, (Dorodnitsyn et al. 2008 b). 



This is clearly seen in the O VIII La line which is observed in absorption at almost all 
inclinations at which warm absorber is seen, and all times. Also, virial arguments suggest 
that the absorber in this case is located closer to the torus. At later times, when the 
torus has the shape of a thick disk, at higher inclinations we are still able to see gas close 
to the source, where Vp ~ Vip, and the M-shaped profile is observed. Unfortunately, the 
fine structure of most lines is seen only in high resolution spectra and completely washed 
out by the convolution procedure. At most times, three dimensional modeling of our 2.5 
hydrodynamical models predicts warm absorber spectrum in a relatively narrow range of 
inclinations: < 10°. 



4-0.3. Fe line region 

We have calculated spectra from a torus wind in the 6.4 keV region. The results for the 
A3 model are shown in Figure |8] as a function of inclination, 6. Since the fluorescent FeK a 
line cannot be produced by resonance scattering, we do not have it in our spectra. 

There are two prominent absorption features from resonant absorption in FeXXV 1.85A 
(6701 eV) and FeXXVI1.78A (6960 eV). Both lines have v{Xo) ^ 100 kms-^ Blueshifts 
are consistent with that expected from the flow becoming slower at higher inclinations. 
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The centroid energy is located at moderate blueshifts: for FeXXVI, a.t 6 = 35°, f(Ao) 
170kms~^, and v{X'^) = 688kms~^, becoming at = 45°, f (Aq) ^ lOOkms"^, and ^^(A^) 
688 km s-^ 
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Fig. 8. — Model A3, observed at t = 4 a different inclinations. Notice prominent absorption 
lines from FeXXV1.85A (6.701 keV) and FeXXVI1.78A (6.96 keV). 

Analogously, from the FeXXV line we have: at ^ = 35°, f(Ao) = 179kms~^, and 
f(A^) = 671kms~^, becoming at 6' = 45°, v{Xq) = 134kms~^, and v{X~^) = 573kms~^. 

The results for the O VII La line and for Mg XI lines are shown in Figure |9] Thus, from 
O VIII La line we obtain A"^^ = IQA (652.9 eV), A^^^ = 18.95A (654.4 eV) where A"^^/^;^ 
is the energy at the maximum/minimum of the corresponding emission component. From 
this we derive v{\^^^) = — 312kms~^ and v{\^^^) = 376kms~^. Analogously, we obtain 
v{X-) = -725 km s"^ and f(A+) = 789 km s^^ 

Performing the same analysis, for example, for the MgXI line we obtain: v{\:^^^) ~ 
Okms~^ and v{X^^^) = 443kms~^, and for the maximum value of the velocity we obtain 
f(A~) = — 665kms~^ and f (A^) = 887kms~^. Thus it appears for this model the Fe La are 
moving lightly more slowly than the O VIII or Mg XI lines. 
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Fig. 9. — Model A3, observed at t = 4 a. 9 = 45°. Profiles of individual lines of O VII (left) 
and MgXI (right). Blueshifts correspond to positive velocities, redshifts to negative 
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Fig. 10. — Model B^, observed at t = 4 a = 55°. Unconvolvcd spectrum is shown. Notice, 
blending of many lines due to O and low ionization Fe, near 0.8 keV. 
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4-0. 4- Spectrum at high inclinations 



At inclinations, 6' > 50 — 60°, continuum radiation of the nucleus is heavily obscured. 



and only radiation scattered in lines is seen. Figure 10 shows the unconvolved spectrum 
from B3 model at t=4, 6 = 55° revealing numerous emission features. Comparing with 



the prototypical NGC 1068, Syfert 2 spectrum of Kinkhabwala et al. (2002), we can see 



that many lines observed in emission in AGN spectrum can be identified in our simulations. 
However in general, the method adopted in this paper tends to overestimate the radiation 
flux scattered in lines. 



The apparent emission bump seen in the Figure Pi)^ at ~ 16 - 16A is significantly 



more prominent than in observations of Kinkhabwala et al. (2002). Individual lines which 



contribute to this hump such as lines from Ne IX, Fe XIX, Fe XIII, O III, and O VII and 
moderately ionized iron, Fe XVII - Fe XXIV are observed in NGC 1068 and are easily 
identifiable in our simulations. Notice that in our method we treat all lines as resonant. 
This simplification does not lead to any significant errors when calculating the absorption 
spectrum, although in emission more radiation is scattered in lines as it would be if the source 



function (12) would include "thermal" term. In the present studies we thus overestimate 



the scattered flux, and a good example of it is the prevalence of the emission about 15A 
compared to real AGNs. 

At inclinations of 70 — 80° the spectra of B models become completely black. As for A 
models, for example, at inclinations as high as 80° at t=4.5 the ^3 has a spectrum, which is 



very similar to that shown in the Figure 10 Overall emission is on the 10% level of that of 



the continuum. Notice an emission excess near 15. 5A (0.8 keV), which we have discussed in 
the context of properties of spectra from A and B models. The spectrum contains numerous 
features from O, Mg, Si, Ne, and low ionization Fe. The existence of such spectrum suggests 
a continuous distribution of ionization parameter, ^ over the scattering region in qualitative 



correspondence with observations of Seyfert 2 galaxies (Kinkhabwala et al. 2002 Brinkman 
et al.||2002l). 



5. Conclusions 



We have carried out 3D radiation transfer simulations of the absorption-emission spectra 
produced by the X-ray excited wind from the obscuring torus in AGN. As a previous stage 
of our research we performed time-dependent 2.5D numerical simulations of wind which 
originates from cold, rotationally supported torus whose inner throat is exposed to a strong 
ionizing radiation of the BH accretion disk (Dorodnitsyn et al. 2008 a|[b ). Previous results 
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include purely absorption spectra from such a flow, and suggest that this model is promising 
in establishing cold tori as a major source of the warm absorber flows observed in many 
AGN. 

Motivated by the fact that in case of an extended wind a ID transmission approximation 
for the radiation transfer in lines is inadequate, we relax this approximation by performing 
3D simulations of the radiation transfer in X-ray lines in a moving gas of the torus wind. 
The distribution of the gas is obtained from the hydrodynamical modeling which is used as 
an input to the radiation transfer calculations carried out making use of a 3D generalization 



of the Castor escape probability formalism (Rybicki & Hummer 1978, 1983) 



The richest warm absorber spectra are observed within ~ 10° — 15°, and the in- 
ferred velocity of the gas decreases as inclination increases. Synthetic spectra predicted by 
our hydrodynamical models are generally consistent with observations of warm absorbers. 
They contain many strong absorption as well as emission lines from O, Mg, Si, Ne, and Fe. 
Numerous features found in our synthetic spectra near 13.7-13.9A (892 - 905 eV), 9.3-9.4A 
(1319 - 1333 eV), and 6.7 — 6. 9A (1797 - 1831 eV) are likely due to absorption from high ion- 
ization Ne, Mg and Si, respectively. Many lines found in our synthetic spectra are blended, 
in particular with those from Fe XVII — Fe XXIV. Velocities of this medium ionization com- 
ponent which are inferred directly from spectral lines shift are consistent with those derived 
directly from hydrodynamical models. Some of the individual line profiles demonstrate P- 
Cygni profiles and others have weak M-shaped profiles, with the red-ward wing being usually 
stronger than the blue-ward wing. Such a profile is the result of complicated flow geometry 
at smaller radii: i.e. a combination of P-Cygni profile and the importance of the azimuthal, 
and poloidal, vg velocities at smaller radii. 

If a warm absorber flow is thermally driven or ablated off of a cold or nearly-neutral 
surface, then the observed low ionization component provides evidence for this origin. This is 
confirmed by our numerical simulations: at the angles at which the line of site goes through 
the X-ray evaporated "skin" of the torus we see such a stripped gas. Observations also suggest 
the simultaneous presence of gas at a low ionization state, log — and a higher ionization 
component with log ^ ~ 100. Note also some claims of observational evidence ruling out 



the existence of gas at intermediate ionization parameters (McKernan et al. 2007; Holczer 



et al. 2007) for some objects. Gas at ionization parameters outside the range 0<log(^) <3 



is harder to detect unambiguously from observations. 

Torus winds naturally possesses features which are necessary for the formation of the 
flow with characteristics typical for warm absorbers. Our calculations show that the angular 
distribution of ^ is strongly angle-dependent, reflecting the 2D nature of the flow. 
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Our solutions also show that ^ can vary dramatically along the line-of-sight. In ab- 
sorption, at such inclinations at which the warm absorber is observed, the line-of-sight first 
touches the stripped cooler gas and then penetrates the torus funnel, where hotter gas is 
located. Emission line photons are collected from a larger volume, coming from separated 
places in the outflowing gas, those which are not significantly obscured by the cold gas of the 
torus. In such flow different photons in a single line may come from regions quite separated 
in space (multiple resonances) making a single zone approximation for ^ inadequate. 

Our numerical methods are not capable of capturing some physics which is probably of 
importance here, most notably the co-existance of multiple phases in the outflowing plasma 
on level smaller than a typical size of our numerical cell. Bearing in mind all these limitations, 
we note that in our hydrodynamic calculations we don't see specific strong bi-modality in 
the distribution of ^. 

Our models divide in two broad categories: more dense, having initial Compton optical 
depth of the torus, = few -10, and less dense, with ~ 1. Additionally, we consider 
cases in which we put the initial torus, at either Ipc, or closer, 0.5pc from the BH. The 
X-ray heated gas is accelerated by the thermal mechanism within the throat of the torus. 
As already established ( Dorodnitsyn et al.||2008 b ) such a funnel mechanism helps the gas to 



accelerate to velocities, Vmax of the order of several escape velocities, V^sc(lpc). This velocity, 
Knax is not necessarily equal to the terminal velocity, V°°, due to the properties of the funnel 
flow. 

Results from our new method generally agree with the pure transmission calculations of 



Dorodnitsyn et al. (2008 b). That is, warm absorbers are observed roughly in a 10° range. 
In general, the velocity corresponding to the centroid energy of an absorption line, is smaller 
than 300kms~^. Most of such absorption is produced in a hot X-ray evaporated skin of a 
torus. At such distances, the finite size of an X-ray source is not important. However, the 
blue-ward part of the absorption line in many cases moves as fast as 1000 km s~^, meaning 
the absorber is at least partially located at smaller radius. At such small radii the core 
can no longer be treated as a point source. Rather it subtends a finite solid angle which 
finally translates to a broader angular pattern at which the absorption spectrum is seen from 
infinity. Such an effect is revealed from individual line profiles. The most rapidly moving gas 
{v > 500 — 800 kms~^) sees such core from close proximity, and that influences the maximum 
blue- and redshift of photons within the line profile. 

This suggests that the torus wind may extend inward to the region where considered 
for accretion disk winds. As it could be impossible to tell between the torus wind and an 
accretion disk wind at such small distances from a BH, this result implies an accretion disk 
wind should be considered as a possibly important ingredient of modeling of warm absorber 
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flows. 
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